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Abstract. Supercurrents in superconducting flux threaded loops are expected to oscillate 
with the magnetic flux with a period of hc/2e. This is indeed true for s-wave superconductors 
larger than the coherence length ifo- Here we show that for superconductors with gap nodes, 
there is no such strict condition for the supercurrent to be hc/2e rather than hc/e periodic. 
For nodal superconductors, the flux induced Doppler shift of the near nodal states leads to 
a flux dependent occupation probability of quasi-particles circulating clockwise and counter 
clockwise around the loop, which leads to an hc/e periodic component of the supercurrent, 
even at zero temperature. We analyze this phenomenon on a cylinder in an approximative 
analytic approach and also numerically within the framework of the BCS theory. Specifically 
for rf-wave pairing, we show that the hc/e periodic current component decreases with the 
inverse radius of the loop and investigate its temperature dependence. 



PACS numbers: 74.20.Fg, 74.25.Fy, 74.25.Sv 
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1. Introduction 

Electrons moving in a multiply connected geometry threaded by a magnetic flux d> are an 
ideal system to observe quantum mechanical phase coherence. If an electron encircles a flux 
threaded hole on a closed path, the phase difference of its wave function must be a multiple 
of 2n plus the Aharonov-Bohm phase 2n<£>/<&o where <£>o = hc/e is the flux quantum [1]. 
Therefore, a finite phase gradient persists in the wave function for all flux values 0/<5o + 1 ■ 
Consequently, a persistent current is flowing around the hole, which is modulated by the 
magnetic flux with a period of <J> [2, 3]. 

This phenomenon is best observed in phase coherent superconducting (SC) rings [4, 5, 6]. 
Measurements of magnetic flux trapped in SC rings showed that the flux is quantized in 
multiples of Oo/2 [7, 8], which implies a flux periodicity of hc/2e for the circulating 
supercurrent and likewise for all thermodynamic quantities [9]. Indeed, the same periodicity 
has been found by Little and Parks in measurements of the critical temperature T c of flux 
threaded cylinders [10, 11]. 

The hc/2e flux periodicity of the SC state is naturally contained in the BCS pairing 
theory of superconductivity [12], as was shown by Byers and Yang [5] and independently by 
Brenig [13] and by Onsager [14]. Byers and Yang introduced two distinct classes of SC wave 
functions, which are not related by a gauge transformation; one class of states has minima in 
the free energy at even multiples of the SC flux quantum fl>o/2, whereas the second class has 
minima at odd multiples of <J>o/2. They proved that the minima in the free energy become 
degenerate in the thermodynamic limit and all thermodynamic quantities hc/2e periodic. 
In finite systems, this degeneracy is lifted. Consequences for flux-dependent oscillations 
of T c have been investigated by Bogachek et al [15] using a quasi one dimensional (ID) 
thin-ring model with s-wave pairing. This and other recent works, including an analysis of 
the supercurrent, made evident that s-wave rings smaller than the SC coherence length 
display a flux periodicity of hc/e rather than the anticipated hc/2e periodicity [16, 17, 18, 19]. 
Special attention arose from numerical investigations of SC c/-wave loops [20, 21], in which 
the distinction of the two classes of SC states is more pronounced than in s-wave rings. These 
results agree with the analytical approaches using the thin-ring model [22, 23] and will be 
extended in this article towards a two dimensional (2D) multi-channel model. 

For the discussion of the flux periodicity of the supercurrent we choose a discrete 2D 
lattice in a cylindrical geometry (figure 1). We recall that the magnetic flux threading a 
superconducting loop is quantized in units of the superconducting flux quantum hc/2e [7, 8]. 
This quantization reflects the minima of the free energy [5]. These minima are determined 
by gauge invariance and the electron interaction; the flux quantum hc/2e is therefore a 
fundamental property of any superconductor. Flux quantization in a cylinder requires that 
its walls are thicker than the penetration depth A. If the walls are thinner than A, the cylinder 
can be threaded by an arbitrary magnetic flux and only the quantity called fiuxoid is quantized 
[4, 9, 6]. In this situation, it is the flux periodicity of thermodynamical quantities such as the 
supercurrent or T c , for which the pairing of electrons suggests hc/2e periodicity. 

For two reasons we expect nodal rather than nodeless superconductors to support an 
hc/e periodicity. The first arises from the discrete nature of the eigenenergies in a finite 
system. The results of the summation over occupied eigenstates for integer and half-integer 
flux values differ by an amount proportional to the mean level spacing Sf in the vicinity of 
the Fermi energy E F . In the normal state, 6 F oc l/V, where V is the volume of the system. 
For the thin cylinder shown in figure 1 with a circumference Na and a height Ma, where a 
is the lattice constant, the level spacing is 8f °c 1/(NM); in s-wave superconductors with 
an order parameter A » Sf, 6f matters little. For SC states with gap nodes, the situation is 
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Figure 1. As a model system to study persistent supercurrents we use a thin-wall cylinder 
constructed of a 2D discrete lattice. The interior of the cylinder is threaded by a magnetic flux 
<I>; we assume that the flux does not penetrate into the cylinder itself. In such a system, <D can 
be chosen arbitrarily, since quantization applies to the fluxoid and not the flux itself. 

different. For example, in the of-wave superconductors with an order parameter Ak oc - k\, 
the nodal states closest to Ef have to fulfill the condition k z = k^, thus there are fewer possible 
eigenstates and Sf cc 1 /N. 

The second reason is that for gapless superconductors with a finite density of states 
(DOS) close to E F , the occupation probabilities of these states change with flux. The flux 
dependence of the occupation enhances the difference of current matrix elements for integer 
and half-integer flux values [20, 19, 23]. This effect can be understood in terms of the spacial 
extension of a Cooper pair. In s-wave superconductors, the occupation probability remains 
constant for all i>, if the diameter of the cylinder is larger than £o- If this condition is fulfilled, 
the constituents of a Cooper pair cannot circulate separately; the pair does not "feel" the 
multiply connected geometry of the cylinder. But for nodal SC states, the length scale which 
characterizes their coherence, diverges in the nodal directions and there are always Cooper 
pairs which extend around the circumference of the cylinder. Therefore nodal superconductors 
have no characteristic length scale above which the SC state is unaffected by the geometry of 
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the system. These two combined effects are investigated on the basis of an analytical model 
in section 3 and by numerical calculations in section 4. 

2. Superconductivity in a Flux-Threaded Cylinder 

The properties of a finite-size multiply connected superconductor depend sensitively on the 
discrete energy spectrum in the normal state, in particular in circular symmetric geometries. 
To understand the SC spectrum of the discrete NxM lattice, we therefore have to characterize 
first its normal state spectrum. To illustrate the problem, we consider the tight-binding 
spectrum of a ID ring with N lattice sites and nearest-neighbor hopping t. A half-filled band 
corresponds to chemical potential p — in equation (2). If N/4 is an integer, there is an energy 
level at energy e = for (f> = 0, where (p = 0/<[>o is the dimensionless magnetic flux. If N/4 is 
a half-integer, the levels are symmetrically distributed above and below e = (figure 2). As a 
function of flux, the spectrum is /z/e-periodic in both cases. If N is odd, there are two possible 
configurations of energy levels, as shown in figure 2(c). In both configurations with odd N, 
two levels cross E F in one flux period. The combination of a particle-like and a hole-like band, 
used to construct the SC spectrum, then becomes hc/2e periodic. These number dependent, 
qualitative differences control the flux dependence of the normal persistent current, as was 
shown by Biittiker et al and by Cheung et al [2, 24]. 

Whenever an energy level crosses Ef with increasing flux, the current reverses its sign, 
thus it is /z/e-periodic for even N and either paramagnetic or diamagnetic in the vicinity 
of (p = 0, and it is /zc/2e-periodic for odd N. The lattice-size dependence persists also in 
rings with electron-electron interactions [25, 26, 27] or in mesoscopic SC islands [28] and in 
particular in a 2D cylinder geometry with circumference Na and height Ma. Each energy level 
of the ID case splits up into M levels, which results in a characteristic flux dependence of the 
spectral density. For special ratios N/M, the flux values where the ID levels cross have a high 
degeneracy; for N - M, the degree of degeneracy is M. For the latter case, the differences 
between the spectrum for integer and half-integer flux values are most pronounced; they are 
similar to the ID spectrum of figure 2(a), if N and M are even, and similar to the spectrum 
of figure 2(b), if N and M are odd. For N = M ± 1, the spectrum is almost /zc/2e-periodic, 
which is the extension of the odd N case in the ID ring. Away from these special choices 
of N and M, the degeneracies are lifted, indicated by the blue shaded "patches" in figure 2. 
The inclusion of a next-nearest neighbor hopping term or a change of p in equation (2) has 
a similar effect, as shown by Zhu [21]. The size of the normal persistent current circulating 
around the cylinder is controlled by the change of the DOS near Ep upon increasing 0. Since 
normal persistent currents in metallic rings are typically hc/e periodic [1, 29], we will choose 
N - M and p = for our model study, where the hc/e periodicity of the spectrum is most 
clearly established, and we will use even N and M for all subsequent calculations. In this 
section and in section 3 we show how this size dependent features survive into a SC state with 
gap nodes. 

The starting point for our investigations is the BSC-theory formulated on a flux threaded 
cylinder with circumference Na = 2nRa and height Ma, where R is the dimensionless radius 
of the cylinder and a the lattice constant. The pairing Hamiltonian is given by 



where k = (k^,,k z ) with k v = n/R and n e {-N/2 + 1, . . . ,N/2). In the z-direction along 
the axis of the cylinder, we choose open boundary conditions, which allow for even-parity 
solutions with k z = (2m e - \)n/M and odd-parity solutions with k z = 2nm /M, where 




(1) 
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Figure 2. The energy spectrum of a cylinder in the normal state depends on the numbers N and 
M, which parametrize the circumference and height of the cylinder. The black lines represent 
the energy levels for a ID ring with M = 1 and (a) N/4 an integer, (b) N/4 a half integer and 
(c) N an odd number. In (a) and (b), level crossings occur at each multiple of the maximum 
Doppler shift for <p = 1/2, denoted by li (see section 3). For odd N, two different spectra 
without level crossings at E(</>) = are possible [N = 4n + 1 (left) and N = 4n - 1 (right) with 
n £ N]. At the red points, a level crosses the Fermi energy Ef = 0. For M » 1, the levels split 
up and form a quasi continuous DOS that depends on the ratio N/M (blue patches). 



m e ,m e {1, . . . ,M/2}. The operators ct and c k . s are creation and annihilation operators for 
electrons with crystal angular momentum hn and crystal momentum hk z Ja. The eigenenergies 
of free electrons moving on a discrete lattice on the surface of the flux threaded cylinder have 
the form 



e k (0) = -It 



*(**-!) 



cos L + cos k- 



H- (2) 



For R » 1, ek(0) can be expanded to linear order in <p/R and 

%(0-%(O)*-2t|sinV (3) 
K 

is commonly called the Doppler shift. 

The superconducting order parameter in the pairing Hamiltonian (1) is defined through 

A k (q, <p) = A q ((p)g(k) = i Y V(k, k')<c k<T c_ k - +qi - c kl c_ k , +qT > (4) 
L k' 

where V(k, k') is the pairing interaction. Here we choose the interaction in the separable form 
V(k, k') = Vg(k)g{k') with the pairing interaction strength V. The order parameter A k (q, </>) 
represents spin-singlet Cooper pairs with total crystal angular momentum hq. On the cylinder, 
the coherent motion of the Cooper pairs is possible only in the azimuthal direction, therefore 
q = (q/R,0) with q e {-N/2 + I,... ,N/2}. The quantum number q is chosen to minimize 
the free energy. The ^-dependence of A f/ (0) enters through the self-consistency condition and 



Flux Periodicities in Loops of Nodal Superconductors 



6 




Figure 3. Scheme for the evolution of the multiply degenerate energy levels in the normal 
state with increasing order parameter A for the rf-wave SC state. For <f> = and q = (a), there 
is a degenerate level E(k, 0, 0) = that splits up for finite A into levels spreading between -A 
and A. For 0=1/2 and q = 1 (b), there is an energy gap around E = of width 2l\, which 
persists into the SC state. 



has been discussed extensively in [16] and [19] for s-wave pairing, where g(k) = const. Since 
A ? (0) varies only little with <f>, we shall postpone the discussion of the flux-dependence of 
the li-wave order parameter to the numerical evaluations of section 4 and start our analytical 
calculation with a <p and q independent order parameter A CJ (<p) = A. As in our preceding 
work [19], we take q = floor(2</> + 1/2) in a first step, such that <p - q/2 is hc/2e periodic; 
eventual deviations from this relation will be discussed in section 4. Since the Hamiltonian 
( 1 ) is invariant under the simultaneous transformation <f> — > (p + 1 and q — > q + 2, it is sufficient 
to consider q = or 1 and the corresponding flux sectors -1/4 < <f> < 1/4 and 1/4 < < 3/4, 
respectively. 

The diagonalization of the Hamiltonian (1) leads to the quasi-particle dispersion 

£ ± (k, q, 4>) = — ± yjAl + e 2 (k, q, <f>), (5) 

with e(k, q, </>) = [ek(</>) + e-k+ q (</>)]/2. Expanding £ ± (k, q, </>) to linear order in both <p/R and 
q/R gives 



£ ± (k,q,0) * -e,(k) ± ^A2+(e k (0)-/ 9 (k)) 2 



(6) 



where 



— q/2 tq 

e 9 (k) = 2tsink (f and / ? (k) = — sin^. (7) 

R R 

In the normal state A = 0, the additive combination of e q (k) and l q (V) leads to the q- 

independent spectrum equation (2). For A > 0, the spectrum (6) differs for even and odd q, 

except for special ratios of and M, as discussed above. This difference is crucial for nodal 

SC states, as shown schematically in figure 3 (and especially for af-wave pairing in figure 9): 

The condition ss k z for levels close to £> causes a level spacing 6f ~ 2/i (kjr) for small A, 

where kf is the Fermi momentum. For and M even and q = 0, the degenerate energy level 

at E — Ef = splits into M levels for increasing A, which spread between -A and A. For 

q = 1, the degenerate levels closest to Ef are located at E = ±|/i(kf)|, thus a gap of 2/i(k/r) 

remains in the SC spectrum. If N and M are odd, the spectra for even and odd q (figure 3(a) 
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and (b)) are exchanged, and if either N or M is odd, the spectrum is a superposition of (a) and 
(b). 

The gauge invariant circulating supercurrent is given by 

Jq(.<f>) = \ X v k«k. s (q), (8) 

where vt = dek(<f>) / d{Rk,p) is the group velocity of the single-particle state with eigenenergy 
fk(<A)- The spin independent occupation probability of this state is 

«b(q) = <cL c k. s >(q) = « 2 (k, q, <f>)f(E + (K q, 0)) - v 2 (k, q, 0)/(E_(k, q, 0)), (9) 

where /(£) = 1/(1 + e E!kBT ) is the Fermi distribution function for the temperature T. The 
Bogoliubov amplitudes are 



o 1 
t/(k, q, 0) = 2 



e(k,q,0) +1 



and v 2 (k,q,0) = 1 - M 2 (k,q,</>). (10) 



From equations (8) and (9), the supercurrent in the cylinder is obtained by evaluating 
either the sum numerically, as discussed in section 4, or from the approximative analytic 
solution in section 3, which allows insight into the origin of the /zc/e-periodicity in nodal 
superconductors. 



3. Analytic Solution and Qualitative Discussion 

An analytic evaluation of the supercurrent is possible only in the thermodynamic limit where 
the sum over discrete eigenstates is replaced by an integral. For a multiply connected 
geometry, this limit is not properly defined because the supercurrent or the Doppler shift 
vanish in the limit R — > oo. Care is needed to modify the limiting procedure in a suitable way 
to access the limit of a large but non-infinite radius of the cylinder. In this limit it is mandatory 
to consider the supercurrent density j q (<f>) = Jq(<f>)/M rather than the supercurrent Jq(<p). In 
this scheme, we treat the DOS as a continuous function in any energy range where the level 
spacing is oc l/NM, but we keep the finite energy gap of width 2/ ? (k F ) oc I /Roc l/N around 
E F in the odd-q sectors. 

For a tight binding energy spectrum as defined in equation (2), the density of states is 
a complete elliptic integral of the first kind. For the purpose of an analytic calculation, a 
quadratic spectrum with a constant density of states in two dimensions is a more appropriate 
starting point. We use the expanded form of equation (2): 



e k (<f>) = t 



Mr 



(ID 



where//' = p+4-t. The quadratic spectrum (1 1) has no upper bound and the sum in equation (8) 
correspondingly extends from -oo to oo for both k v and k z . 

Some algebraic steps are needed to rearrange the sum in equation (8) suitably to convert it 
into an integral. For finite <p, ey^cfi) + £-k(<P), and consequently the sum has to be decomposed 
into a component with > and a second one with < 0. We therefore take k^ > and 
write Vk as 

v ±k = | [+K - |) = v d (k) ± v,(k), (12) 

where v^(k) = -2t(p/R 2 is the diamagnetic contribution and v p (k) = Itk^/R, a paramagnetic 
contribution, respectively [30]. 
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In a continuous energy integration, the Doppler shift is noticeable only in the vicinity of 
E F . On the Fermi surface, and k z are related by: 



In this spirit we approximate e q (k) and l q (k) by e q (k z ) « 2r(0 - q/2)k^ F (k z )/R and Z 9 (fc z ) « 
tqk v jr(k z )/R, respectively. The eigenenergies (6) near £> are thereby rewritten as 

E + (±k^, k z , q, cf>) = +e q {k z ) + ^Al + (e k (0) + wf ( 14 ) 

E_{±k v , k z , q, 0) = - ^A 2 k + (e k (0) + l q (L)f (15) 

For the evaluation of the supercurrent J q ((/>) in equation (8), the sum over k is now 
replaced by an integral over k^ and k z , which is then performed by integrating over the normal 
state energy e and an angular variable 9. According to our scheme for replacement of discrete 
energy levels by a continuous spectrum, the DOS becomes gapless in the limit M — > oo for 
q = 0, although N is kept finite. For q — 1 instead, a ^-dependent gap 2|/i(fc z )| remains. Thus 
we replace 6^(0) + \l q (k z )\ by the continuous quantity e± \l q (E F , 9)\. In summary, the procedure 
is defined by the following steps: 

) — > — dk v dk z = — dkk d6 = MN de dO, (16) 

V 2n -JO 27T Jo J-n/2 J-u' J -Ti/2 



~ JO JO J-n/2 J -ft 1 J-n/2 

where we use the parametrization 



k v \ k cos 9 \ / e + // I cos 9 



k z J \ ksinO J \ t \ sin 9 



(17) 



with e = tk 2 - p.' and where N - R/Ant is the constant DOS in the normal state. The energy 
integral runs over the whole tight-binding band width 8f with the Fermi energy E F = in the 
center of the band. Correspondingly, we integrate from -//' to //. Furthermore, the Doppler 
shift is parametrized for e « E F as 

e q (9) = ^~^ /2 2tylJ7Ttcos6 and Z,(0) = ^ -^Tfcos 0, (18) 

where the function / 9 (0) is positive for all allowed values of 9. The supercurrent thus becomes 

1 e r 



J ** )= Mh 



^ Vk«k s (q) + ^ v k n ks (q) 

* 2A/| J | d0 J" de[n ?+ (e, 0)v + (e, 0) + n q -(e, 0)v_(e, 9)1 (19) 

where n q± {e,6) = n ± k(e,e)(q) an d v ± (e,8) = v ±k(e,e)- The factor 2 in equation (19) originates 
from the spin sum. We collect the terms proportional to v d (e, 9) = -2t(f>/R 2 into a 
diamagnetic current contribution jj and the terms proportional to v p (e, 9) = 2tk^ F (e, 9)/R 
into a paramagnetic contribution j p . Using f(-E) = 1 - f(E), equation (19) simplifies to 

•Vr/2 

(20) 



j" = 4N if de r de ^ ^ / 1 \f (E + - k- e + - 

fcJ-*/2 J W fl) VA 2 + e 2 L J 
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where ja = jd(q, 4>) and j p = j p (q, <f>). Here, the integration is over positive values of e only 
and the lower boundaries of the integration over e are controlled by l q (ff). Since l q (k) = 
at the minimum of the band (e = -//'), the upper integral boundary remains p' . We used 
the abbreviations A = A(6) and E = E(e, 6) = ^A 2 {9) + e 2 . The current j d turns out to be 
diamagnetic in the even-g flux sectors and paramagnetic in the odd-q sectors. For even q, it 
is equivalent to the diamagnetic current obtained from the London equations [31, 32]. The 
current j p has always the inverse sign of j d and is related to the quasi-particle current as shown 
below. As presented in section 2, E displays a distinct spectra in the even-g and odd-q flux 
sectors. To analyze the flux dependent properties of the spectra and the current, we distinguish 
the case of s-wave pairing (or any other superconducting state with a complete energy gap) 
and the case of unconventional pairing with nodes in the gap function. For the latter, we focus 
on of- wave pairing. 



3.1. s-Wave Pairing Symmetry 

For s-wave pairing, A(e, 6) = A is constant. Therefore, if we assume that A > e q (9) for all 6, 
the lower energy integration boundaries in equations (20) and (21) are equal to A. Thus j q (<p) 
is equal in both the even-g and odd-q flux sectors and the flux periodicity is hc/2e. However, 
if A < max e e q (9), equation (8) has to be evaluated exactly, the procedure and results of which 
have been presented in [19]. 

With e = V£ 2 - A 2 , equations (20) and (21) transform into integrals over E with 
de = D S (E) dE, where 

Dm _ ^ _ J E(E 2 -A 2 r m for E>A 

Ds(E) ~ 8E - I for £ < A (22) 



is the SC density of states for s-wave pairing. This leads to 

-vr/2 



jd = 4N lf_ de £ dEv " ( V£ 2 - A 2 , 9) [f(E + e q {8)) - f(-E + e q (6))] , (23) 
j P =W- h [' d0 f dED s (E)v p ( VZ? 2 - A 2 , 0) [f(-E - e q {9)) - f(-E + e q (0))] . (24) 



nil 

At T = 0, we find 



— d9 cos 6 dED s {E)+0[-\ . (26) 

t J -nil J A V 1 1 



StNe 
~R~h 



>-n/2 

\2 



In the integral over of j p , the inequality e/t <sc 1 applies, and terms of order 0(e/t) 2 are 
negligible. 

The current jd becomes independent of the SC density of states. Its size is essentially 
proportional to £>, as long as p' » A holds. The paramagnetic current j p depends on the 
absolute value of the order parameter and on its symmetry. 

If A > e q (ff) for all values of 6, then j p = and the supercurrent j q = j^ consists of 
the diamagnetic part only. For T > 0, jd decreases slightly, but remains of the same order 
of magnitude. The current j p increases with increasing T and reaches its maximum value 
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Figure 4. Scheme for the density of states of a d- wave superconductor for <j> = 1/4, where 
£o(l/4) = h/2. The center of mass angular angular momentum hq of the Cooper pairs is 
(a) q = and (b) q = 1 . The energies are Doppler shifted to higher (red) or lower energies 
(blue). This results in a double-peak structure and for q = in an overlap of the upper and 
lower "band" in the region -\e$\ < E < \eo\ [34] and states in the upper band become partially 
occupied. For q = 1 there is a gap / 1 of the size of the maximum Doppler shift at <j> = 1/4. The 
black line represents the density of states (a) for <f> = and (b) for <f> = 1/2. 



at T c . For finite temperatures j p is usually denoted as the quasi-particle current. The entire 
supercurrent is always the sum of the diamagnetic current and the quasi-particle current 
j p , and therefore decreases with temperature and vanishes at T c [33]. The quasi-particle 
current has the same flux periodicity as the supercurrent, even though it is carried by single 
quasiparticle excitations. In the normal state (A = 0), 



Jp 



8tNe 
~R~h 



, «r/2 pW) ,ed>-q/2 r' 2 

— d6cos6 dE = 4p'-- — — , 

t J-n/2 JO h Rn J_ n/2 



dQ cos 2 6 



•9/2 



h R 

which cancels jj exactly in the limit M 



(27) 



3.2. Unconventional Pairing with Gap Nodes 

For a more general order parameter A(0), an analytic solution of equations (20) and (21) is 
hard to obtain. For s-wave symmetry, depends only weakly A; is indeed maximal for 
A = 0. Equation (25) for j c { is valid also for unconventional order parameter symmetries. 
Physically, jd reflects the difference in the DOS of quasi-particle states with orbital magnetic 
moment parallel and anti-parallel to the external magnetic field. The first group of states is 
Doppler shifted to lower energies, whereas the latter is Doppler shifted to higher energies. 
This is schematically shown in figure 4 for c/-wave pairing (c.f. [34]). In this picture, 
jd is proportional to the difference between the area beneath the red and and blue curves 
representing the DOS arising of bands £L(±|k|,q, <p) < (underlaid red and blue). Therefore 
we approximate jd for A(ff) <c p! — E p + At by 



% In this procedure, the normal persistent current vanishes, but this is unproblematic here because the normal current 
above T c is exponentially small for k%T c » r5f . 
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as given by equation (25) with A = 0. On the other hand, j p is represented by the occupied 
quasi-particle states in the overlap region of E + (k, q, </>) and £_(k, q, </>) with width 2e q (k F ). It 
is therefore strongly dependent on the characteristic DOS in the vicinity of E F . In figure 4(a), 
which refers to even q, the current j p is determined by the small triangular patch where the 
upper and lower bands overlap. For odd q, the two bands do not overlap, therefore j p = 0. 

We will now analyze such a scenario for d-wave pairing. With an order parameter 
A k = A(k 2 - k 2 ) a A cos 26. Again, we assume A > e q (ff) for all 6; then the integral in 
equation (21) contains only the nodal states closest to E F , for which the d-wave symmetry 
demands « k z . Jointly with equation (13) this condition fixes the Doppler shift at Ep to the 
k-independent value e q = (0 - q/2) y/2tfi'/R and l q = (q/R) y/tfi' /2. With the density of states 
in the d-wave superconducting state 

D d (E) = — I (29) 
V£ 2 - A 2 cos 2 28 

equation (21) for the paramagnetic current j p at T = then takes the form 
e r« r 12 2t je + u' 

jp=4N h{ ^Ij^ni-r^ 6 - (30) 

In the odd-<7 flux sectors, l q > e q for all values of 0, therefore j p = 0. In the q = sector, 

2e I/I' r« ,„ r 12 ,„ . n i 2e fi7 r« E e fi7 2 

nhA V t q 



l q = and 



2e hi' C e « , f" z , . 1 2e hi' f c " , E 

jp~—\— dE dOsmG dE- = 

hn V t Jo J-x/2 V£ 2 - A 2 cos 2 29 nh \ t J A 



2 L~73 e (<P-i/2 



2 



nA^h\R)' (31) 
where the same approximations as in the s-wave case are applied. The dominant contribution 
to the angular integral ocer 6 originates from the nodal parts, where the integrand can be 
linearized in 6, such that the integral can be performed approximately (see e.g. [35]). 
In the even-g sectors, the total current j q = jj + j p finally becomes 

V^ 7 ' 



1 



7TA R 



(32) 



which results in the ratio 

j p V^ 7 4> 



= b<f> (33) 



id r 

of the two current components. 

In the odd-g flux sectors, j p = and the supercurrent is j q (0) = jd- As a function of (/>, 
j q ((/>) is consequently hc/e periodic; within one flux period from -1/2 to 1/2 we represent it 

as 

, [0+1/2 for -1/2 < <p < -1/4, 
M) = -2^ e -{ 0(1 - H) for -1/4 < <f> < 1/4, (34) 
K 11 { <p - 1 12 for 1/4 < < 1/2, 

(c.f. figure 5). The amount by which the supercurrent differs in the even-g and odd-g flux 
sectors is represented best in the form of Fourier components: the n-th Fourier component of 
j((p) is j„ = d(f> j(<p)e 2mn<1 ' . Here, we denote the first Fourier component by jh/ e , and the 
second Fourier component by jh/2e and obtain 

u! e ,8-tt 2 . u'eAni-b 

]hle = - 2 Rh b ^ ^ Jhl2e = - 2 R-h^- (35) 
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Figure 5. The supercurrent density j q (<p) = jd + j p in a thin t/-wave cylinder as a function 
of flux <f> (arbitrary units). Shown is the result of the analytical model (equation (34)) for the 
characteristic value b = 0.4. For -1/4 < <p < 1/4, where q = 0, the current is reduced by a 
contribution proportional to ifr, whereas it is linear in (p otherwise. This gives rise to an overall 
flux periodicity of hc/e. 



To leading order in l/R, the ratio of the h/e and the h/2e Fourier component therefore is 



Jh/e 



Jh/2e 



0.07- 



2t 
AR 



(36) 



and scales with the inverse ring diameter. This 1 /R law is the direct consequence of the d- 
wave density of states Dd(E) cc E. For some other unconventional superconducting states with 
D(E) cc E" in the vicinity of Ef, the decay of the jh/ e Fourier component results in a 1 /R" 
law. Using equation (36) to estimate this ratio for a mesoscopic cylinder with a circumference 
Ra = 2600a ~ 1 ^m and a ratio A/t = 0.01, we obtain ji,/ e / jh/2e ~ 0.03. 



4. Numerical Solution for ti-Wave Pairing at T = 

In this section we evaluate numerically the supercurrent in equation (8) together with the 
self-consistency condition 

£ = X . 8kiq) [/(£-*', q, 0) - 4>)] , 07) 

where the tZ-wave pairing symmetry follows from by 

g k (q) = cos(^ - q/2) - cos k z (38) 

and the order parameter is Ak(q, 0) = A q ((p)g^(q). Here we take into account the full q- and 
^-dependence of Ak(q, <f>). The g-dependence of gkiq) is essential to ensure the invariance of 
the gap equation (37) under the replacement (f> — > <p + 1 an d 9 ~ > 1 + 2. At those flux values 
for which the total energies 

E t (q,(f>) = Y j e k ((f>)n ks (q). (39) 

are equal, q advances to the next integer. This flux value may deviate from the values 
= (2n - l)/4, for which we fixed the g-sector transitions in section 2. 

Loops of c/-wave superconductors can be arranged in two different ways. In a first choice 
for the geometry the order parameter windes jointly with the lattice around a hole such that the 
phase of the order parameter remains constant on the selected path. The cylinder geometry 
described here is an example for this choice. The second option is to fix the orientation of 
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Figure 6. Energy spectrum and supercurrent J q ((f>) in a cylinder with circumference of N = 56 
and height M = 28 and Ao = \(4> = 0) as 0.1?. The red lines represent occupied states 
which form the condensate, whereas the black lines represent empty states. The spectrum (a) 
is similar to the one obtained in [20] for a square frame geometry. Clearly visible is the energy 
gap in the odA-q flux sectors, whereas in the even-g flux sectors states cross the Fermi energy 
upon changing the flux. At these crossing points, a jump in the supercurrent is observed (b). 



the lattice and to cut out a hole. Then the phase of the order parameter rotates by 2n on any 
closed path encircling the hole once. This was investigated with a square frame in [20] and 
also in the ID model in [22]. These two arrangements are in fact physically equivalent. The 
square frame geometry ensures the right number of lattice sites for the maximum difference 
in the spectrum of the even-q and oAA-q flux sectors, as discussed in section 2. For a direct 
comparison to the cylinder geometry, we chose a cylinder with N - 56 and M = 28, which 
has the same hole diameter as the square frame in [20], and the ratio N/M = 2 produces 
qualitatively the same energy spectrum. The resulting spectrum is shown in figure 6(a). It has 
indeed the same characteristic features as in the square frame geometry. An energy gap of 
the same order of magnitude exists in the odd-q flux sectors, and the DOS in the even-g flux 
sectors is gapless. There are no hybridization effects in the spectrum of the cylinder, since it 
preserves the full rotational symmetry. The features mentioned above are also in agreement 
with the qualitative discussion of section 3. Clearly visible in figure 6(b) are the jumps in the 
supercurrent whenever an energy level crosses Ef, and the offset in the flux value for which 
q changes (large jumps). This offset depends in a complex way on the system size and the 
pairing potential strength, but generally decreases for larger values of N, M, and V. 

The spectrum and the supercurrent in figure 6 display the expected signatures of 
discreteness which are not captured by the analytic analysis of section 3. The important 
parameter is obviously the size of the level spacing. Explicitly we take a closer look at 
a cylinder with N - M - 2600, and thus a circumference of the order of 1 /an. The 
calculated spectra are shown in figure 7(a) and (b) for different pairing potentials, resulting 
in (a) A q (<p) * 0.05f and (b) A q (<p) * 0.02f. The qualitative features ot the much smaller 
cylinder remain, but the gap /; in the odd-q flux sectors is smaller, because /; decreases with 
1 IN. In the even-g flux sectors, there are M levels spread out between -A q ((p) and A q (<f>), 
which leads to an increase in the DOS around Ef with decreasing A q (<f>) for fixed N and M. 
The representation with a continuous DOS is therefore appropriate for AJ(p) <K t, which is 
fulfilled well in figure 7(b). Figure 7(c) shows the spectrum for M - 1 = N = 2600, which 
is almost identical in even-g and odd-q flux sectors. There is still a gap for non-integer (or 
half-integer) values of <f>, but it is equally distributed in the even-g and odd-q sectors. Other 
choices of N and M produce mixed features of the spectra in figure 7(a) and 7(c). All the 
energy levels shown in each part of figure 7, which belong to nodal states, have apparently the 
same derivative with respect to 0. 
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(c) V = OAt 

Figure 7. A section of the energy spectrum around E F . (a) and (b) show the spectra of a 
cylinder with N = M = 2600. (a) V = OAt and A (0) * 0.05f; (b) V = 0.2t and A (0) * 0.02f. 
The energy gap l[ <K A 9 (0) for these systems and all the states shown have the same Doppler 
shift (all lines are parallel / perpendicular). The density of states is quasi continuous in the 
even-q flux sectors and grows linearly with decreasing A q (<f>). (c) shows the Ac/2e-periodic 
spectrum of a cylinder with N = 2600 and M = N + 1. 

The small level spacing in the pm sized cylinders results in solutions A q (<p) of the gap 
equation (37), which are nearly constant (figure 8(a), note the vertical scale discussed in the 
figure caption). The ^-dependence of the total energy E,(q, (p) also becomes small, whereas 
the small difference for even-g and o&A-q remains important for the supercurrent J q ((f>). Since 
J q (</>) oc dE,(q, <p)/d</>, the differences in E,(q, (p) imply different current amplitudes in the even 
and odd q sectors [see figure 8(b) and (c)]. This effect is larger for smaller A q {<f>), because the 
number of energy levels crossing E F increases with decreasing A q ((f>). For the chosen pairing 
potentials V, the difference of the amplitudes of J q ((f>) for even and odd q are of the order of 
a few percent. Per contra, the current jumps within a ^-sector are tiny for the large radius of 
the /Lmi-size cylinders. However, the resulting A q {<p) is considerably larger than in the d-wave 
cuprate superconductors. § Consequently, the upper limit for the difference of J q ((p) in the 
even-<7 and odd-g sectors would be larger in the cuprate superconductors than in the model 
system calculated here. The offset in the jump of J q ((/>) from the flux values (p = (2n - l)/4 
is resolved for the smallest V in figure 8. However, at low temperatures the superconducting 
state for each q becomes meta stable for those flux values, for which it is not the ground state. 
At which flux values such a meta stable state decays into the ground state is not clear and the 

§ Angle-resolved photoemission spectroscopy on various cuprates suggests a tight-binding t x 200 meV - 400 meV. 
The gap at the antinodes, obtained from tunneling spectra, varies between 10 meV and 50 meV [36, 37], therefore 
A » 0.002* -0.01?. 
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Figure 8. (a-c) Self-consistent order parameter A g (0), total energy E t (<j>) and supercurrent 
J q (f) plotted as a function of flux <p for V = 0.25? (dark blue), V = 0.2? (light blue), and 
V = 0.15? (green). A q (<f>) and E,(<p) are shown in units of ? and J q (<j>) in units of et/hc. The 
oscillations of all quantities are oc l/R and of the order of ? for E,(<p). The amplitude of the 



oscillations in A q (<p) are rather small. For V 
amplitude <SA = [A (0) - A,(l/2)]/A (0) » 10"* 



: 0.25?, A q (cj>) x 0.036? with an oscillation 
A q (<f>) a 0.02? and SA*5x 10~ 6 (V = 0.2?), 
and A q (<t>) a 0.009? and <SA ss 4 X 10~ 5 (V = 0.15?). (d) Ratio of the first and second Fourier 
components of the supercurrent as a function of the cylinder radius R for fixed values of A. The 
height M of the cylinder is equal to N = 2nR, which yields the maximum values for Ji,/ e . For 
N larger than some A-dependent number (see main text), the results of the exact evaluations fit 
very well to the prediction of equation (36) (black line). 



position of the jump in the supercurrent can vary in experiments. 

We now compare the ^-dependence of the ratio of the first and second Fourier 
components Jh/e/Jh/2e analogous to section 3. This is shown in figure 8(d) for different values 
of V. The ratio is in excellent agreement with equation (36) for system sizes larger than a 
few hundred lattice constants. For smaller systems, Jh/ e becomes larger than predicted by the 
1 /R size dependence. The scale which decides about the validity of the approximations used 
in section 3 is the ratio of the level spacing and A q (ip). Equation (36) therefore holds, if the 
prefactor of <p in equation (33) is small, that is, if V8 t «: nAR, because t/R is proportional 
to the level spacing of the nodal states. For a cylinder with radius Ra = 2600a « 1 pm and 
Aq{4>)/t ~ 0.01, we obtain the ratio Jh/e/Jh/2e ~ 0.04, which is almost identical with the result 
of section 3. 

5. Periodicity Crossover for small A 

So far we always assumed that A 9 (0) » e q {(j)) and concluded that variations in A q {<p) are 
negligible. But if A q {<p) is of the same order as the Doppler shift e q {<p), the situation changes 
dramatically. This is the case, if either the radius R of the cylinder is very small, or the pairing 
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Figure 9. Periodicity crossover at T = in a cylinder with N = M = 400 for V = 0.40? (dark 
blue), V = 0.21 (light blue), V = 0.15/ (green), V = O.lf (orange) and V = O.Olt (purple). For 
V = 0. 1 1 and V = 0.07f, A g (0) is smaller than the Doppler shift for all <f>, thus J q ($) approaches 
the /ic/e-periodic normal persistent current. 



potential V is small or the temperature T is close to T c . Here, we analyze the flux periodicity 
and the crossover from a "small-gap" to a "large gap" regime by increasing V from zero to 
higher values at T = 0, and by lowering T through T c for fixed V. As mentioned above, 
the amplitude of the oscillations, especially those of A q {(f>), become very small for increasing 
R. For very large R, the periodicity crossover takes place within a tiny range of V or T, 
respectively. To observe the crossover more comfortably in a larger window of V or T, we 
use smaller systems here. 

The mechanism of the periodicity crossover at T — 0, controlled by V, is best discussed 
by analyzing the total energy E,(q, <f>) [figure 9(b)]. It differs little from the crossover in s-wave 
superconductors, for which it was investigated in [19]. In the normal state (V = 0), E,(q, <p) is 
^-independent and consists of an hc/e periodic series of parabolae. For increasing V, a new 
minimum in E,(q, <p) forms at the crossing points of two parabolae. This minimum mover 
downward in energy until this new parabolic arc crosses the neighboring parabolae at the flux 
values <p = (2n - l)/4. The energies of the old and the new minima are generally different for 
any finite system, but they approach each other when A q {<p) » 6f . In the odd-q flux sectors, 
A q {(p) is nearly constant because no energy levels cross E F , whereas in the even-g sectors, 
levels cross £> for all values of V. This causes the wiggles in E,(q, <p) and the decrease of 
A q (tp) with increasing (p [figure 9(a)]. For the smallest two values shown in figure 9(a), A ? (0) 
approaches zero as a function of <p for even q; for this reason the odd-g states extend far into 
the even-g flux sectors. With increasing V, the nearly hc/2e periodic sawtooth pattern of the 
supercurrent evolves from the hc/e periodic normal persistent current [figure 9(c)]. 

The temperature controlled crossover at T c is analogous to the crossover controlled by V, 
but the finite temperature has quenched all the effects of discreteness as well as the gap in the 
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Figure 10. Temperature driven periodicity crossover for fixed V = 0.4? in a cylinder with 
N = M = 100 for k B T = 0.1863? (dark blue), k B T = 0.1870? (light blue) and k B T = 0.1873? 
(green). The amplitude of the normal persistent current in the sectors with A q (<p) = is much 
smaller than for A c/ (0) > and is invisible on this plot scale. 



odd-<7 flux sectors. This means that the deviations from the hc/2e periodicity are invisible 
in figure 10. Deviations appear with decreasing temperature as k B T approaches l\. The 
supercurrent decreases linearly with increasing T until it reaches the exponentially small value 
of the normal persistent current at T = T c [figure 10(b)] [33]. This suppression as well as the 
suppression of A q (<p) with temperature [figure 10(a)] differ only little from those of s-wave 
superconductors. The only qualitative difference is, that a characteristic temperature T* exists 
for s-wave superconductors, below which A{<p = 0) is larger than the maximum Doppler shift. 
This is equivalent to a coherence length f(T*) = 2R [19]. Below T", A(<p) > for all <p in 
i-wave superconductors, and the thermodynamic quantities are therefore not affected by the 
Doppler shift. The relation A(0 = 0, T — 0) ~ 1.75T C leads to the estimate 

T -^- « Ef (40) 

For af-wave pairing, there is no such characteristic temperature because of the nodal states, 
but in analogy we can define T* as the crossover temperature below which A 9 (0) > for all 
cf>. Analogously to the s-wave case, we denote this situation as the "large-gap" regime. For 
temperatures T* < T < T c , A 9 (0) approaches zero for certain values of </>, which we call 
the "small-gap" regime. Since for a of-wave superconductor with nearest neighbor hopping 
A(<p = 0, T = 0) > 1.75£b?c [38], one expects that T c - T* is also larger and the crossover 
broader than for s-wave pairing. 



6. Conclusions 



We have shown that in rings of unconventional superconductors with gap nodes, there is a 
paramagnetic, quasi-particle-like contribution j p > to the supercurrent at T — 0. This 
current is generated by the flux induced "reoccupation" of nodal quasiparticle states slightly 
below and above E F . Formally a coherence length hvf/A^iq, <p) > 2R can be ascribed to 
these reoccupied states, which are therefore affected by the geometry of the system, however 
large the number of lattice sites is. If the normal state energy spectrum has a flux periodicity 
of hc/e, than the superconducting spectrum is hc/e periodic, too. The normal state spectrum 
of a cylinder with a discrete lattice strongly depends on the number of lattice sites on the 
cylinder. This problem is characteristic for rotationally symmetric systems and is much less 
pronounced in geometries with lower symmetry, such as the square frame discussed in [20]. 
In such systems, the addition or removal of a small number of lattice sites or impurities do not 
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change the spectrum qualitatively, as tested by numerical calculations on a square frame. For 
an experimental arrangement where the difference in even and odd flux values is as large as 
possible, a square loop would be preferable. Our results obtained in section 3 and section 4 
for the periodicity of the physical quantities A 9 (0), E,(q, (/>), and J q ((f>) provide therefore an 
upper limit for the hc/e periodic components. 

The hc/e periodicity is best visible in the current component j p at T = 0. For d- 
wave-pairing j p oc 1/R 2 , and the hc/e periodic Fourier component decays like the inverse 
radius of the cylinder, relative to the hc/2e periodic Fourier component. The lack of a 
characteristic length scale in nodal superconductors, such as the coherence length for s-wave 
pairing, generates this algebraic decay with increasing R. Although j p is larger for small A, it 
almost vanishes close to T c , if A » 6f , and variations of T c with flux, as in the Little-Parks 
experiment [10, 11], do not differ for s- and <f-wave superconductors. 

A possible set-up for the experimental detection of the hc/e periodicity of the 
supercurrent is the insertion of Josephson junctions onto the cylinder, thereby creating a 
SQUID. The oscillations of the SQUID's critical current have the same flux periodicity as 
the circulating supercurrent. Indeed, experiments with cZ-wave SQUIDs by Schneider and 
Mannhart have shown an hc/e periodic Fourier component under certain conditions [39]. The 
relation to the effect described here however is not established, yet because of the so far 
unexplored influence of the Josephson junctions. 

A different approach to study the cross over from the normal persistent current to the 
supercurrent in a ring was proposed by Biittiker and Klapwijk [40] and later by Cayssol et al 
[41]. They analyzed a normal metal ring with an s-wave SC segment of variable length I. The 
energy spectrum, which they found, depends on I in a similar way as it does in our analysis 
on the radius R. In this set-up, hc/e periodicity should be found if I < £ > although the ring 
diameter is much larger than £o- Analogously, we expect the ratio j\/ e / jh/2e to be proportional 
to 1 // for a J-wave SC segment. This might be of advantage for experimental detection. 
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